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Abstract 

The theoretical description of the correlations between observables in two separated rapidity 
intervals in relativistic nuclear collisions is presented. It is shown, that the event-by-event pt~Pt 
correlation defined as the correlation between event mean values of transverse momenta of all particles 
emitted in two different rapidity intervals does not decrease to zero with the increase of the number 
of strings in contrast with two particle pt~Pt correlation - the correlation between the transverse 
momenta of single particles produced in these two rapidity windows. 

In the idealized case with the homogeneous string distribution in the transverse plane in the 
framework of the cellular analog of string fusion model (SFM) the asymptotic of pt-pt correlation 
coefficient is analytically calculated and compared with the results of the Monte-Carlo (MC) calcula- 
tions fulfilled both in the framework of the original SFM and in the framework of its cellular analog, 
which enables to control the MC algorithms. 

In the case with the realistic nucleon distribution density of colliding nuclei the results of the MC 
calculations of the pt-pt correlation function for minimum bias nuclear collisions at SPS, RHIC and 
LHC energies are presented and analysed. 

1 String fusion model (SFM) 

The colour string model [1, 2] originating from Gribov-Regge approach is being widely applied for the 
description of the soft part of the multiparticle production in hadronic and nuclear interactions at high 
energies. In this model at first stage of hadronic interaction the formation of the extended objects - 
the quark-gluon strings - takes place. At second stage the hadronization of these strings produces the 
observed hadrons. In the original version the strings evolve independently and the observed spectra are 
just the sum of individual string spectra. However in the case of nuclear collision, with growing energy 
and atomic number of colliding nuclei, the number of strings grows and one has to take into account the 
interaction between them. 

One of possible approaches to the problem is the colour string fusion model [3] . The model is based 
on a simple observation that due to final transverse dimensions of strings they inevitably have to start 
to overlap with the rise of their density in transverse plane. At that the interaction of string colour 
fields takes place, which changes the process of their fragmentation into hadrons as compared with the 
fragmentation of independent strings. So we have one more interesting nonlinear phenomenon in nuclear 
interactions at high energies - the field of physics the investigations in which were initiated by pioneer 
works of academician A.M. Baldin [4]. 

It was shown [3, 5, 6] that the string fusion phenomenon considerably damps the charged particle 
multiplicity and simultaneously increase their mean p t value as compared with the case of independent 
strings. In accordance with a general Schwinger idea [7] and the following papers [8, 9] (colour ropes 
model) two possible versions of string fusion mechanism were suggested. 

The first version [5] of the model assumes that the colour fields are summing up only locally in the 
area of overlaps of strings in the transverse plane. So we will refer to this case as a local fusion or overlaps. 
In this case one has 

(n) fe = Mo — Vk (p 2 t ) k =WVk (1) 

Here (n) k is the average multiplicity of charged particles originated from the area Sfe, where k strings 
are overlapping, and {p 2 ) k is the same for their squared transverse momentum. The fig and p 2 are the 
average multiplicity and squared transverse momentum of charged particles produced from a decay of 
one single string, and <7o is its transverse area. 

In the second version [10] of the model one assumes that the colour fields are summing up globally - 
over total area of each cluster in the transverse plane - into one average colour field. This case corresponds 



to the summing of the source colour charges. We will refer to this case as a global fusion or clusters. In 
this case we have 

(n) cl =Vo — Vk7i (p?>c«=pV^ k cl = ^f^ (2) 

Here (n) cl is the average multiplicity of charged particles originated from the cluster of the area S c i and 
(Pt) cl is the same for their squared transverse momentum. The N^ r is the number of strings forming 
the cluster. 

2 Cellular analog of SFM 

To simplify calculations in the case of real nucleus-nucleus collisions a simple cellular model originating 
from the string fusion model was proposed [11]. In the framework of the cellular analog along with the 
calculation simplifications the asymptotics of correlation coefficients at large and small string densities 
can be found analytically in the idealized case with the homogeneous string distribution, which enables 
to use these asymptotics later for the control of the Monte-Carlo (MC) algorithms (see below). 

Two versions of the cellular model as the original SFM can be formulated - with local and global 
string fusion. In this model we divide all transverse (impact parameter) plane into sells of order of the 
transverse string size Co- 

In the version with local fusion the assumption of the model is that if the number of strings belonging 
to the ij-th cell is kij , then they form higher colour string, which emits in average fio \fk~Tj particles with 
mean p\ equal to p 1 '\fk~ij (compare with (1)). Note that zero "occupation numbers" kij = are also 
admitted. 

In the version with global fusion at first we define the neighbour cells as the cells with a common link. 
Then we define the cluster as the set of neighbour cells with non zero occupation numbers kij ^ 0. After 
that we can apply the same formulae of the global fusion (2) as in the original SFM, where N^ r is the 
number of strings in the cluster and S c i/a n is the number of cells in the cluster. 

From event to event the number of strings kij in the ij-th cell will fluctuate around some average 
value - k^ . Clear that in the case of real nuclear collisions these average values kij will be different for 
different cells. They will depend on the position (sjj) of the ij-th cell in the impact parameter plane (s 
is two dimensional vector in the transverse plane). In the case of nucleus-nucleus AB-collision at some 
fixed value of impact parameter b one can find this average local density of primary strings kij in the 
point Sij using nuclear profile functions Ta{s^ + b/2) and Ts{sij — b/2). 

In MC approach knowing the kij one can generate some configuration C = {kij}. To get the physical 
answer for one given event (configuration C) we have to sum the contributions from different cells in 
accordance with local or global algorithm (see above), which corresponds to the integration over s in 
transverse plane. Then we have to sum over events (over different configurations C). Note that as 
the event-by-event fluctuations of the impact parameter at a level of a few fermi are inevitable in the 
experiment one has to include the impact parameter b into definition of configuration C = {b, kij}. 

3 Long-range correlations 

The idea [5, 6, 12] to use the study of long-range correlations in nuclear collisions for observation of the 
colour string fusion phenomenon based on the consideration that the quark-gluon string is an extended 
object which fragmentation gives the contribution to wide rapidity range. This can be an origin of the 
long-range correlations in rapidity space between observables in two different and separated rapidity 
intervals. Usually in an experiment they choose these two separated rapidity intervals in different hemi- 
spheres of the emission of secondary particles one in the forward and another in the backward in the center 
mass system. So sometimes these long-range rapidity correlations are referred as the forward-backward 
correlations (FBC). 

In principle one can study three types of such long-range correlations: 
n-n - the correlation between multiplicities of charged particles in these rapidity intervals, 
Pt-pt - the correlation between transverse momenta in these intervals and 
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p t -n - the correlation between the transverse momentum in one rapidity interval and the multiplicity of 
charged particles in another interval. 

Usually to describe these correlations numerically one studies the average value (B) F of one dynamical 
variable B in the backward rapidity window Ay B , as a function of another dynamical variable F in the 
forward rapidity window Ay F . Here (...} F denotes averaging over events having a fixed value of the 
variable F in the forward rapidity window. The (...) denotes averaging over all events. So we find the 
correlation function (B) F = f(F). 

It's naturally then to define the correlation coefficient as the response of (B) F on the variations of 
the variable F in the vicinity of its average value (F) . At that useful also to go to the relative variables, 
i.e. to measure a deviation of F from its average value (F) in units of (F), and the same for B. So 
it's reasonable to define a correlation coefficient b B _ F for correlation between observables B and F in 
backward and forward rapidity windows in the following way: 



(F) d(B) F 



bB ~ F ~ (B) dF 



(3) 

F={F) 



As the dynamical variables we use the multiplicity of charged particles (n), produced in the given 
event in the given rapidity window, and the event (!) mean value of their transverse momentum (p t ), i.e. 
the sum of the transverse momentum magnitudes of all charged particles, produced in the given event in 
the given rapidity window (Aj/), divided by the number of these particles (n): 

1 ™ 

Pt= ~^2\Pti\, where y { e Ay; i = l,...,n. (4) 

i=i 



So we can define three correlation coefficients: b n _ n , b pt _ pt and b pt _ n , for example: 

(n F ) d(p tB ), lF 

i u p t -n — 

P tF = {p tF ) 



(p tF ) d( PtB ) Ptf 
~ ( PtB ) dp tF 



^ ~ (PtB> dn F 



(5) 

™ F =<™ F > 



The n B , n F are the multiplicities and p tB , p tF are the event (4) mean transverse momentum of the 
charged particles, produced in the given event correspondingly in the backward (Ay B ) and forward 
(Ay F ) rapidity windows. In this paper we concentrate only on the pt~Pt - correlation. Some other results 
can be found in [13]. 



4 Analytical results and MC calculations 

In the framework of the cellular analog of SFM in the idealized case with the homogeneous string dis- 
tribution in the transverse plane (kij = r\ = const for all i and j) the asymptotic of pt-pt correlation 
coefficient can be found analytically at large string density p str = v/ao, n >• 1, when one cluster is 
forming: 

, _ fJ-OF Ah- VOF ( fi\ 

pt - pt Mof + 16 7 2 Vl ~ vl ^f + 47 2 VV { ' 

Here /jqf is the average number of charged particles produced from a decay of one string in the forward 
rapidity window. The 7 is the coefficient of proportionality between the average transverse momentum 
p and the square root from the dispersion of p for one string: <r p = jp, where a p = p 2 —p 2 . The 7 is 
the same for p 2 : a p 2 = ^yp 2 . In the last formula of (6) as the dynamical variables in the forward (F) and 
backward (B) rapidity windows instead of p tF and p tB (4) we have used: 

n F ^ n B 

PtF = P ^F ' aild = P ^B . ( 7 ) 

U F i= i n B i=1 



The asymptotics (6) are similar to the ones obtained in [11] for n-n and p t -n correlations: 
u u Vof i7 1 Mof 

6 — = ^-" = ^T^ and 6 -" = 2^T4^' (8) 
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notwithstanding the fact that the derivation of (6) is very different and much more complicated. We see 
that for pt-pt correlation in asymptotic the /xo/y^-scaling obtained in [11] for n-n and p t -n correlations 
is also takes place. The correlation coefficients depends only on the one combination of /xo and r\. 



Cellular analog of SFM 
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Figure 1: The comparison of the asymptotic of pt-pt correlation coefficient b pt _ pt (6) with the results 
of its MC calculations in the framework of the cellular analog of SFM in the idealized case with the 
homogeneous string distribution in the transverse plane at two values of /iyF = A*o • The rj is proportional 
to the string density 77 = (JoP str - The total number of sells is M = 450. See the text for details. 
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Figure 2: The same as in Fig.l but for MC calculations in the framework of the original SFM at the 
average number of strings N str =1000, 4000 and 8000. 

In Fig.l we compare the asymptotic (6) of b pt _ pt with the results of its MC calculations in the 
framework of the cellular analog of SFM described above. We see that at r\ > 5 the results for both 
local and global versions of the model practically coincide with the asymptotic. The MC calculations 
in Fig.l fulfilled for the total number of cells M — 450, but the simulations with the small number of 
cells M = 45 give practically the same results. So for pt-pt correlation one also has "M-scaling", which 
was found earlier in [11] for p t -n and n-n correlations. This means that the pt-pt correlation coefficient 
depends only on the string density p — n/o~o but not on the total number of strings N str = nM. 

In Fig. 2 we compare the same asymptotic (6) with the results of MC calculations in the framework 
of the original SFM (with the taking into account the real geometry of merging strings, which takes 
much more efforts). Again we see that at w > 5 the MC results calculated for local (overlaps) and global 
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(clusters) versions of the original SFM practically coincide with the asymptotic (6). In Fig. 2 the points 
correspond to the MC calculations fulfilled at the different number of strings N str =1000, 4000 and 8000. 
We see that in this case the results also do not depend on the total number of strings but depend only 
on the string density p str = rj/ao- 

We would like to emphasize that this takes place only due to the fact that as observables p tF and p tB 
we choose the event mean values of transverse momenta of the particles emitted in the rapidity intervals 
(see (4) and (7)). In paper [14] it was shown that if one takes as observables the transverse momenta 
of single particles produced in these two rapidity windows (two particle correlation) then such defined 
Pt~Pt correlation coefficient is inversely proportional to the number of strings N str and rapidly decreases 
to zero with the increase of the total number of strings. 

Comparing the Fig.l and Fig. 2 we see also that the MC calculations in the framework of the original 
SFM and in the framework of its cellular analog give very similar results at all values of string density. 

5 The pt-pt correlation in minimum bias AA-collisions 

In the case with the realistic nucleon distribution density of colliding nuclei the results of the MC cal- 
culations of the pt~Pt correlation function (see section 3): (p tB ) p = f(p t F) f° r minimum bias nuclear 
collisions at SPS, RHIC and LHC energies arc presented in Fig. 3. Certainly, in the case of minimum bias 
MC calculations one has to include the impact parameter b into definition of configuration C (see section 
2). In these figures (o) and (•) denote the results of calculations in the framework of the original SFM 
for its local (overlaps) and global (clusters) versions correspondingly. The (□) and (■) denote the results 
of calculations in the framework of the cellular analog of SFM for its local and global (clusters) versions. 
All presented results are for the forward rapidity window of 2 unit length (Ay F = 2). 

First of all we see in Fig. 3 that the results of MC calculations for all four versions are practically 
coincide with each other. Only the results for both versions with local fusion arc slightly shifted to the 
higher values of p t in comparison with the ones for clusters, which can be easily explained by simple 
general arguments. 

We see also that in the case of minimum bias nuclear collisions the pt~Pt correlation functions (p tB ) p 
exhibit strongly nonlinear behavior. Note that the correlation functions calculated at fixed value of impact 
parameter (or with its small fluctuations within a certain centrality class) are practically linear and the 
strength of the correlation can be described by one number - the correlation coefficient (see (5)). It's not 
true for minimum bias collisions and the correlation can't be described by one number in this case - the 
derivative of the correlation function is positive at small values of p tF and it becomes negative at large 
values of p tF . 

This can be explained in the following way. The set of minimum bias collisions includes both central 
and peripheral collisions. In our model with string fusion due to higher density of strings one has higher 
value of p tB and p tF for central collisions, than for peripheral ones. This is the reason for the positive 
Pt~Pt correlation. To understand the reason for the negative pt~Pt correlation at large values of p tF one 
has to remember that as observables p tF and p tB (4) (7) we choose the event mean values of transverse 
momenta of the particles emitted in the rapidity intervals. Only in this case the pt~Pt correlation is not 
equal to zero in real heavy ion collisions (see the discursion in section 4). It's clear from (4) (7) that 
although the event mean values of p tF and p tB are higher for central collisions, their dispersions is higher 
for peripheral ones, as the number of charged particles n F and n B produced in these rapidity intervals 
are smaller for peripheral collisions. This means that if in some event one observes a very high (or a very 
low) value of p tF in the forward rapidity window (than the typical value (p tF )), then it's more probably 
that the peripheral collision has taken place and then one finds the p tB in the backward rapidity window 
at a low level - typical for peripheral collisions. So with the increase of p tF the function {p tB ) p p starts 
and finishes at a low (typical for peripheral collisions) value of p tB obtaining in the middle more higher 
(typical for central collisions) value. 

Using these arguments and taking into account the increase of string density, typical values of p tF and 
p tB , and the multiplicities n F and n B , with the increase of the initial energy, one enables, also to explain 
the change of the form of the pt~Pt minimum bias correlation function in Fig. 3 from SPS to RHIC and 
LHC energies. 
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6 Conclusion 



The theoretical description of the correlations between observables in two separated rapidity intervals in 
relativistic nuclear collisions is presented. The calculations are fulfilled in the framework of the model 
taking into account the possibility of quark-gluon string fusion and also in the framework of the suggested 
simple cellular analog of this model. For both models two possible mechanisms of a string merging are 
considered: the local fusion ("overlaps") and the global fusion (with forming of colour "clusters"). 

It is shown, that for the pt~Pt correlation the event-by-event correlation between event mean values 
of transverse momenta of the particles emitted in two different rapidity intervals does not decrease to 
zero with the increase of the number of strings in contrast with the correlation between the transverse 
momenta of single particles produced in these two rapidity windows which was studied earlier [14]. 

In the idealized case with the homogeneous string distribution in the transverse plane the asymptotic 
ofpt-pt correlation coefficient analytically calculated in the framework of global string fusion mechanisms, 
when at large string density one cluster is forming. The asymptotic is compared with the results of the 
MC calculations in this idealized case both in the framework of the original SFM and in the framework 
of its cellular analog, which enables to control the MC algorithms. 

In the case with the realistic nucleon distribution density of colliding nuclei the results of the MC 
calculations of the pt~Pt correlation function for minimum bias nuclear collisions at SPS, RHIC and LHC 
energies are presented and analysed. 
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Figure 3: The results of the MC calculations of the pt~Pt correlation function (p tB ) p = /(f? t F) f° r 
minimum bias nuclear collisions at SPS, RHIC and LHC energies. See the text for details. The results 
are presented in units of pt\, where pt\ = p is the average transverse momentum of particles produced 
from a decay of one string. 
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